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Inelastic neutron scattering measurements are reported for the quantum antiferromagnetic ma- 
terial Cu2(C5Hi2N2)2Cl4 (CuHpCl). The magnetic excitation spectrum forms a band extending 
from 0.9 meV to 1.4 meV. The spectrum contains two modes that disperse throughout the a — c 
plane of the monoclinic unit cell with less dispersion along the unique b-axis. Simple arguments 
based on the measured dispersion relations and the crystal structure show that a spin ladder model 
is inappropriate for describing CuHpCl. Instead, it is proposed that hydrogen bond mediated ex- 
change interactions between the bi-nuclear molecular units yield a three-dimensional interacting spin 
system with a recurrent triangular motif similar to the Shastry-Sutherland Model (SSM). Model- 
independent analysis based on the first moment sum rule shows that at least four distinct spin pairs 
are strongly correlated and that two of these, including the dimer bond of the corresponding SSM, 
are magnetically frustrated. These results show that CuHpCl should be classified as a frustration 
induced three dimensional quantum spin liquid. 



I. INTRODUCTION 

Organometallic magnetslil are excellent model systems 
in which to explore the intricate auantum many-body 
physics of interacting spin systems .0 They are attractive 
because their energy scales are well matched to efficient 
experimental probes of magnetism and because a wide 
range of magnetic phases are found in these materials. In 
addition to supporting frozen jnagnetic states with fer- 
romagnetism, ferX|imagnetism,Q antiferromagnetism and 
meta-magnetism,EJ organometaUic magnets alssHDrovide 
intriguing examples of quantum spin liquids.QB These 
can be defined as strongly correlated states of interact- 
ing spin systems where time reversal symmetry persists 
at temperatures far below the characteristic energy scale 
for interactions. 

Quantum spin liquids are most commonly found in 
quasi-one-dimensional antiferromagnetic systems such 
as the. JHiiform spin-1 chain J^l tiie alternating spin 
chain ,ErO and the spin-ladder .E3"li3 Low dimensional- 
ity generall5j_£!nhances the phase space for low energy 
fluctuationSjEZi and in one dimension, the effect iS|-to pre- 
clude a frozen state involving Heisenberg spinsu How- 
ever, there are also examples of spin-dimer systems 
where the singlet ground state associated with antifer- 
romagnetically interacting spin pairs survives the ef- 
fects of weaJifiijjriter-dimer interactions in two and three 
dimensions .EJiE^ 

GcometricaLfrustration is an alternate route to stmng 
fluctuatiausQ and there aije theoretical predictionsEj of 
quantumnircj and classicalES spin liquids based Oftjtkis ef- 
fect. While_rpaterials that approximate kagom£3i23 and 
pyrochloreE3E2l antiferromagnets generally have a freez- 
ing transition at sufficiently low temperatures, there are 
other more complex structures where geometrical frustra- 
tion stabilizes an isolated singlet ground state in the T = 



limit.lErE3 SrCu2(B03)2 is a three-dimensional version 
of the so-called Shastry-Sutherland modelEj where the 
intcr-dimer interactions are almost as strong as the intra- 
dimer interactions, and yet they fail to induce gapless 
magnetic fluctuations and long range fEdef because no 
such state can satisfy all interactions St£l Frustration 
is also central to stabilizing the quasi-two-dimensional 
cooperativcj— singlet state in the organometallic mag- 
net PHCCEa While in^Cu2(B03)2 correlations exist 
only within spin-dimcrsj£3 the correlated spin clusters in 
PHCC involve at least 8 spins of which two spin pairs 
are frustrated and provide a positive contribution to the 
ground state energy. PHCC was previously identified as 
an alternating spin chain based on susceptibility, mea- 
surements. This and other misidentificationalini3 indi- 
cate that conventional bulk measurements cannot reli- 
ably determine the origin and nature of spin systems 
with a gapped excitation spectrum and that more so- 
phisticated techniques should be applied to explore these 
unique systems. 

One magnet with a gapped excita- 
tion spectrum that has received considerable attention 
is Cu2(C5Hi2N2)2Cl4p(ipii9(l,4-diazacycloheptane)2Cl4 , 
or for short CuHpCl) .ETO Measurements of the magne- 
tization, magnetic susceptibility, and specific heat show 
that this system has a spin gap A « 0.9 meV, a magnetic 
bandwidth of approximately 0.5 meV, and a saturation 
field Hc2 = 13.2 T. Based on these measurements and 
the crystal structure, it was proposed that CuHpCl is 
a two-leg spin ladder composed of coupled dimers, with 
the intra-dimer bonds of strength Ji = 1.14 meV forming 
the rungs, and inter-dimer bonds of strength J2 ~ 0.2Ji 
forming the legs of the ladders .t^ Subsequent experimen- 
tal results have generally been interpreted in terms of 
this model, and the system has_been the inspiration for 
a number of theoretical studies.Ej"E2 
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FIG. 1. The molecular formula unit of CuHpCl featuring 
two copper atoms in 4+1 square pyramidal coordination£3 
The apical Cu-Cl bonds are shown in yellow. The Cu coordi- 
nation pyramids share an apical edge and have parallel basal 
planes (red lines). 



Despite this extensive body of work, the true nature 
of the spin interactions in CuHpCl has never been con- 
clusively demonstrated. Susceptibility and specific heat 
data are equally well described by spin-iadder, alternat- 
ing chain, and coupled-bilayer models.Ej Some evidence 
that the two-leg ladder model might in fact not be appro- 
priate for CuHpCl canjie.from a previous inelastic neutron 
scattering experimented on a powder sample, where we 
found that the wavevector-integrated magnetic scatter- 
ing intensity did not show the characteristic van Hove 
singularities expected for the magnetic density of states 
of a one-dimensional spin system. In this measurement, 
the wavevector dependence of the energy-integrated mag- 
netic scattering was also inconsistent with the predic- 
tions of the ladder model. However, conclusions from 
this experiment were limited by the coarse wavevector 
resolution of the instrument employed, and the large, 
non-magnetic background signal from the hydrogenous 
sample. 

In this paper, we report inelastic neutron scattering 
measurements performed both on a deuterated powder 
with improved wavevector resolution and on hydroge- 
nous single crystals. The powder and single-crystal mea- 
surements are each independently inconsistent with the 
spin-ladder model. The powder data indicate that the 
strongest dimer bond is different from that of the lad- 
der model. The single crystal measurements show the 
presence of two modes within the 0.5 meV bandwidth, 
a feature not predicted by the spin-ladder model. Con- 
sideration of the measured dispersion relation and the 
structure of the system lead to the conclusion that the 
network of significant magnetic interactions in CuHpCl 
is three dimensional. Model-independent analysis of the 
powder and single-crystal data based on the first mo- 
ment sum-rule provide the contributions of each crys- 
tallographically distinct spin pair to the ground state 
energy. Two classes of spin pairs yield positive contri- 
butions, indicating that geometrical frustration plays an 
important role in stabilizing the quantum spin liquid in 
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TABLE I. Cu-Cu Bond lengths and fractional coordinates 
for CuHpCy-iCalculated from previously determined atomic 
coordinatesEj and measured low temperature (T = 4K) lat- 
tice constants. The lattices of interacting spins formed by the 
bonds are shown on Fig. ^ and Fig. 14. The last column shows 
the contribution of each spin pair to the ground state energy. 
Bonds with the same numerical index are closely related in 
terms of bond vectors and chemical environments such that 
their contributions cannot and need not be distinguished in 
this experiment. 



CuHpCl. 



A. Structure of CuHpCI 



CuHpCl is monoclinic with space group P2i/c and 
room temperature lattice constants a = 13.406(3) 
b = 11.454(2) A, c = 12.605(3) A, and (3 = 115.01(2)°.eI 
The low-temperature lattice constants measured with 
neutron scattering at T = 4 K are a — 13.35(4) A, 
6 = 11.24(6) A, c = 12.72(4) A, and /3 = 115.2(2)°. 
Figure |l| shows the approximately centro-symmetric bin- 
uclear molecular unit containing the spin pair denoted 
by bond 1 in Table 1. Each Cu^+ ion is in a (4-1-1) 
square pyramidal coordination with a CI atom at the 
apex and the base formed by two N and two CI atoms. 
Within the Cu(/i-Cl)2Cu complex of CuHpCl, the co- 
ordinating pyramids share apical edges and have par- 
allel basal planes. Susceptibility measurements for a 
series of compounds with this atomic configuration in- 
dicate intra-molecular exchange constants ranging from 
-1.5 meV (ferromagnetic) to 0.9 meV (AFM) with no 
apparent correlation tJiroughout the range of molecular 
structural parametersO 

The intra-molecular exchange interactions in CuHpCl 
do not create an extended lattice. Because the molecular 
units in CuHpCl interact solely through hydrogen bond- 




FIG. 2. The four different types of near neighbor hy- 
drogen bonded molecular pairs in CuHpCl.LJ Numbers indi- 
cate hydrogen bond distances in A. The orientations of the 
molecules correspond to the lattices of interactions shown in 
Fig. 0. Color coding is the same as for Fig. Fl 



ing, the dimensionality and the overall nature of magnetic 
interactions in this system are entirely determined by hy- 
drogen bonding mediated exchange interactions. Owing 
to the slight deviation of the CuHpCl molecules from 
centro-symmetry, molecular pairs come in two flavors in 
the CuHpCl crystal structure, which we denote by a and 
b. Neglecting flavor distinctions, there are four configu- 
rations for molecular pairs in direct contact as shown in 
Fig. Ij. The sublattices generated by each of these pairs 
are shown in Fig. ^. The numbers on Fig. ^ indicate pos- 
sible spin exchange interactions mediated by hydrogen 
bonding and additional information about each of these 
is listed in Table 1 . Each molecule is part of two molec- 
ular pairs of the type shown in Fig. ^(a). The three dif- 
ferent exchange interactions associated with this molec- 
ular pair are denoted 2, 3, and 10 in Fig. | and Table 
1. Bonds 2 and 3 proceed through a Cu-N-H-Cl-Cu path 
with H-Cl distances ranging from 2.3 A to 2.7 A. Bond 
10 has a similar exchange pathway as for bonds 2 and 
3 but must in addition traverse the entire Cu(/z-Cl)2Cu 
complex. The corresponding exchange interaction should 
therefore be significantly weaker than for bonds 2 and 3. 
Each molecule is also part of four pairs of the type shown 
in Fig. ^(b). Exchange pathways involving Cu-N-C-H-Cl- 
Cu mediate two different exchange interactions that we 
denote 4 and 7, with H-Cl bond distances ranging from 
2.6 A to 2.8 A. Bond 7 should be negligible, however, as 
it involves an apical Cu-Cl contact. The same molecular 
pairs also afford an exchange interaction denoted by 8, 
which proceeds through a Cu-Cl-H-C-H-Cl-Cu pathway 
with H-Cl bond distances in the range 2.6-3.0 A. Each 
molecule also is part of four molecular pairs of the type 
shown in Fig. |2|(c). The exchange pathway involves ei- 
ther Cu-Cl-H-C-N-Cu or Cu-N-H-H-C-N-Cu with H-Cl 
and H-H bond lengths ranging from 2.7 A to 3.1 A. The 
corresponding exchange interactions are denoted 5 and 6. 
Finally, each molecule is also part of two molecular pairs 
of the type shown in Fig. ^(d) with displacement vectors 
±b. The exchange path passes through two 1,4 diaza- 
cycloheptane rings so the corresponding interactions are 
likely to be negligible, especially considering the lower 
coordination number for this interaction. 

The molecular pairs in Fig. ||(a) and (d) yield one- 
dimensional lattices extending along the [101] and [010] 
directions respectively. Inter-molecular interactions cor- 
responding to Figs. 2(b) and (c) on the other hand yield 
b — c and a — b spin planes with surface normals a* and 
c* respectively. Unfortunately it is difficult to predict 
the strength of hydrogen mediated exchange interactions, 
which reportedly can range from 0.1 meV to greater than 
10 meV depending on the chemical environment .e3 Previ- 
ous papers on CuHpCl worked under the assumption that 
the inter-molecular bonds depicted in Fig. ^a) are dom- 
inant, leading to the ladder model shown in Fig. ||(a). 
While we cannot provide firm quantitative information 
about the magnitude of inter-molecular exchange interac- 
tions, we shall present evidence that the 8 inter-molecular 
pairs of the type shown in Figs 2(b)-(c) when combined 



are energetically more significant than those shown in 
Fig. |^(a), and that the intra-molccular spin pair (1) is in 
a frustrated configuration. 



II. EXPERIMENTAL TECHNIQUES 
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FIG. 3. The four different types of lattices generated by 
the inter- molecular interactions shown of Figs, y(a-d). The 
molecules in Figs. ^ are oriented like the spin pairs indicated 
by 1 in this figure. Dashed lines indicate interactions that are 
expected to be very weak. Bond numbers refer to Table 1. 
Sub-lattice color coding coincides with that of Fig. 14. 



The powder sample studied consisted of 5.87 grams of 
deuterated CuHpCl. To produce this sample, perdeuter- 
ated 1,4-diazacycloheptane wip,first prepared following a 
previously published method.Ej The dg-dibromopropane 
and N, A^'-dibenzylethylene-c?4-diamine precursors re- 
quired for this synthesis were prepared from com- 
mercially available IjS-dg-propanediol and ethylene-d4- 
diamine, respectively. The CuHpCl powder was ob- 
tained by rapid cooling from 40° C to 0°C of 1:1 molar 
solutions of anhydrous CuCl2 and the deuterated 1,4- 
diazacycloheptane dissolved in the minimum amount of 
deuterated methanol. Prompt gamma neutron activation 
analysis performed on a portion of the sample confirmed 
95.0(1)% deuteration. 

The single crystal measurements were performed on 
a composite sample with a total mass m « 110 mg. 
This sample consisted of four hydrogenous single crys- 
tals, mutually aligned to within 4.5 degrees. These crys- 
tals were obtained by diffusive grcwath from CuCl2 and 
1,4-diazacycloheptane in methanoLcZl We have produced 
crystals as large as 33 mg using this method. 

Inelastic neutron scattering measurements on both the 
powder and single-crystal samples were performed us- 
ing the SPINS cold neutron triple axis spectrometer at 
the National Institute of Standards and Technology in 
Gaithersburg, Maryland. For the powder experiment, 
the horizontal beam coUimation before the sample was 
50'/ki (A~^) - 80'. Scattered neutrons in the energy 
range 2.6 meV < Ef < 3.7 meV were Bragg reflected by 
a flat, 23 cm wide by 15 cm high pyrolytic graphite an- 
alyzer [PG(002)] at a distance of 91 cm from the sample 
position. The analyzer was followed by an 80' radial col- 
limator and a position-sensitive detector. Cooled Be and 
BeO filters were employed before and after the sample, 
respectively. Data were collected by scanning the scatter- 
ing angle 26 in the range 7° to 114° at fixed incident en- 
ergy Ei. Scans at Ei = 4.0, 4.34, and 4.84 probed inelas- 
tic scattering for energy transfer 0.4 meV < hw < 2.09 
meV and momentum transfer 0.2 < Q < 2.35 A~^, 
with average full width at half maximum (FWHM) res- 
olutions Shiv = 0.14 meV and 5Q = 0.014 k-^lB Back- 
grounds due to incoherent scattering from the analyzer 
and low-angle air scattering of the incident neutron beam 
were measured separately. After subtracting these, the 
data were converted to the normalized scattering inten- 
sity I{Q, fiLo) using the measured incoherent elastic scat- 
tering jfcpm the sample following the procedure detailed 
in Ref.y. 

For the single-crystal measurements, the horizontal 
beam collimation before the sample was 50'/fc^ (A-i) - 



80'. A liquid nitrogen cooled BeO filter was placed after 
the sample and data were collected at fixed final energy 
Ef = 3.7 meV while scanning the incident energy in the 
range 3.95 < Ei < 5.45 meV. A horizontally-focusing 
pyrolytic graphite analyzer with horizontal and vertical 
acceptance angles of 5° and 7° respectively was used in 
conjunction with a single cylindrical detector, which sub- 
tended an angle of 4° in the horizontal plane to an area 
element of the analyzer. In this configuration the aver- 
age instrumental energy resolution for the energy transfer 
range of 0.75 to 1.25 meV was Sfiu) « 0.17(1) meV. Repre- 
sentative values of the projected FWHM wavevector reso- 
lution for the constant-Q scans performed throughout the 
measurement are 5Q\\ = 0.081 and 5Qi_ = 0.065 A^^ 
for the components of the wavevector transfer along the 
principal directions of the resolution ellipse at Q = (100) 
and Tiuj = 1 meV. Scans at constant wavevector transfer 
were performed in both the [hQl) and [hkQ) reciprocal 
lattice planes. 



III. RESULTS 
A. Powder Sample Measurements 

Figure ^ shows the normalized scattering intensity 
/(Q, huj) for CuHpCl at T = 0.3 K and T = 30 K. At 
T ~ 0.3 K, there is a band of inelastic scattering in the 
range of energy transfer 0.9 meV <Tiuj < 1.4 meV, con- 
sistent with our^previous measurements on a hydrogenous 
powder sample.EJ At T = 30 K, the intensity in this Q-Tiw 
range is diminished, which confirms that the correspond- 
ing inelastic scattering cross section is magnetic. 

Figure |^ shows the wavevector averaged scattering in- 
tensity 
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This data is a measure of the magnetic density of states. 
One-dimensional magnets have pronounced van Hove sin- 
gularities at the upper and lower bounds of the magnetic 
excitation spectrum. When convolved with the energy 
resolution function, such singularities would pmduce the 
spectrum shown by the dashed line in Fig. plc^ The in- 
consistency between model and data provides a first in- 
dication that CuHpCl is not a one dimensional spin sy|S=. 
tem. A previous experiment led to the same conclusioncll 
and as it was done on a hydrogenous sample while the 
present sample is deutererated, the comparison of these 
experiments indicates that deuteration does not signifi- 
cantly alter magnetism in CuHpCl. Fig. ||(a) also shows 
that inelastic magnetic neutron scattering from CuHpCl 
is strongly peaked near Qq ss 0.6 A~^, with a secondary 
maxmium near 1.3 A-i. The powder averaged scatter- 
ing intensity from a single-spin dimer with spacing d has 
a peak for Q^d = 1.437r.E3 The data thus indicate sin- 
glet formation in CuHpCl between spins separated by 




Q(A-') 

FIG. 4. Powder inelastic neutron scattering intensity 
i{Q, huj) for CuHpCl at (a) T=0.3 K, and (b) T=30 K ob- 
tained by combining measurements at Ei — 4.84, 4.3, and 
4.0 meV. The figure was produced by binning the data in 
bins of size Shu = 0.03 meV and SQ = 0.016 and 
then coarse-grain averaging to set the effective resolution to 
Ahij = 0.2 meV and 6Q = 0.08 A"^ 

d ~ lA3n/Qo — 7.5 A. This result is also inconsistent 
with the ladder model for CuHpCl, as the dominant bond 
in this model is the intra molecular Cu pair (Fig. |l|) whose 
spacing is only di = 3.376 A. 



B. Single Crystal Measurements 

Inelastic scattering measurements were carried out for 
wavevector transfers at the locations in the (hOl) plane 
indicated on Fig. ||, as well as along the line (IfcO) per- 
pendicular to this plane. Figures show the data so 
obtained, while Figs. |l^ and |l^ summarize the corre- 
sponding dispersion relations derived by fitting the con- 
stant wavevector scans to resolution limited peaks. 

When dynamic correlations are dominated by a single 
dimerized spin pair, there is a well tested RPA theory 
that can accountJEoiLjiiany aspects of the magnetic exci- 
tation spectrumE^Lj As we shall show in the following, 
each spin in CuHpCl takes part in several strongly cor- 
related spin pairs, so the RPA theory is not applicable 
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(100) 

FIG. 6. Locations in the (hOl) plane at which sin- 
gle-crystal inelastic neutron scattering data was obtained for 
CuHpCl. The diameters of the points are proportional to the 
measured first moment of the data, and show that the domi- 
nant satisfied magnetic bond is parallel to [101]. 



in its present form. However, inspection of the excita- 
tion spectra and the crystal structure leads to important 
insights concerning the magnetism in CuHpCl. 

We first note that there are two modes in the mag- 
netic excitation spectrum. This is most easily seen at 
Q = (1,0.5,0) in Fig. 0(c), where two resolution-limited 
peaks are clearly visible at huj = 0.88 and 1.2 meV. 
Two modes are also visible at other wavevectors such 
as Q = (1.167,0,0) in Fig. |(b). In addition. Fig. |(a) 



and Figs. |ll(a-c) show broad or asymmetric peaks that 
are well described by a superposition of two resolution- 
limited Gaussian peaks. 

The spin ladder and alternating chain models for CuH- 
pCl corresponding to Fig. ^(a) and ||(d), respectively, 
both have a single degenerate triplet excitation and are 
therefore inconsistent with the observed two modes. This 
is true despite the fact that Cu(/i-Cl)2Cu complexes 
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FIG. 7. Inelastic neutron scattering data for CuHpCl indi- 
cating less than 0.05 meV dispersion along the (IfcO) direction. 
The scan at (1, 0.5, 0) shows the presence of two modes in the 
excitation spectrum. Solid lines are fits to resolution limited 
gaussians. 



in CuHpCl come with two different orientations as the 
inter-molecular interactions in Fig. |2|(a) and ^d) only 
couple molecules with like orientations. However, the 
inter-molecular interactions in Fig. §(b) and ^(c) link 
molecules with different orientations to form lattices with 
two molecules per unit cell [Figs. §(b) and §(c)]. If the 
ground and excited states maintain the full symmetry 
of the paramagnetic molecule and the spectrum is dom- 
inated by resonant modes, then any lattice except those 
in Figs. 3(a) and (d) has two triplet excited states con- 
sistent with the experimental data. This means that the 
lattices shown in Fig. ^(b) and/or ^(c) are essential parts 
of the interacting spin system in CuHpCl. 

This conclusion is reinforced by Figs. ^ and [l^(b), 
which show dispersion along the (/lOl — h) direction 
[equivalent to {hOh)]. With displacement vectors along 
[101] and [010], the molecular pairs in Fig. 1(a) and |(d) 
cannot yield dispersion along this direction in reciprocal 
space. The implication is again that the inter-molecular 
interactions shown in Figsj2(b) and/or 0(c), and hence 
the lattices shown in Figs|^(b) and/or ^(c), are essential 
parts of the cooperative magnetic network in CuHpCl. 
Figures || and |l2|(c) show that there is also dispersion 
along the a* direction. This implies that the interactions 
corresponding to Fig. ^(b) and ||(b) cannot be the only 
relevant inter-molecular interactions in the problem. Fig- 
ures |ll| and ^(a) show that there is dispersion along the 
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FIG. 8. Inelastic neutron scattering data for CuHpCI show- 
ing dispersion along (/lOO), and perpendicular to the lattice in 
Fig. ^(b). Solid lines are fits to resolution limited gaussians. 
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FIG. 9. Inelastic neutron scattering data for CuHpCl in- 
dicating dispersion along the (hOh) direction. Solid lines are 
fits to resolution limited gaussians. 



c* direction. Consequently the interactions correspond- 
ing to Fig. ^(c) and ^(c) also cannot be the only relevant 
inter-molecular interactions. By inference, at least two 
of the four types of inter-molecular interactions in Fig. ^ 
are important to cooperative magnetism in CuHpCl, and 
the interactions in Figs. ||(b) and/or ^(c) must be part 
of the group. 

Figure 14 shows a projection of the lattice of interac- 
tions in CuHpCl along the b direction. The color-coding 
is consistent with Fig. ^ indicating bonds associated with 
different inter- molecular interactions. As the lattices of 
Fig. ^ interpenetrate, a corollary to the above is that 
CuHpCl is a three-dimensional interacting spin system. 
Nonetheless Figs. and |l3| show that there is less than 
0.05 meV dispersion along Q = (IfcO). Based on the dis- 
cussion above, the absence of dispersion along b cannot 
be due to lack of interactions that couple the system mag- 
netically. A likely alternate explanation is that geomet- 
rical frustration localizes the excitations as has 
served in other geometrically frustrated systems .r^^'^^ 



IV. ANALYSIS 

A. Wavevector Dependent Intensity 

The dynamic spin correlation function S{CI,lu) obeys 
sum rules that can be used to draw additional model in- 
dependent conclusions about magnetism in CuHpCl. The 
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FIG. 10. Inelastic neutron scattering data for CuHpCl 
showing dispersion along the {h, 0,1 — h) direction, which is 
perpendicular the spin ladder in Fig. ^(a). Solid lines are fits 
to resolution limited gaussians. 
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FIG. 11. Inelastic neutron scattering data for CuHpCl. 
(a)-(c) show dispersion along the (101) direction and perpen- 
dicular to the lattice in Fig. ^ (c). Solid lines are fits to 
resolution limited gaussians. 

total moment sum rulell provides an important check on 
whether the measured scattering intensity accounts for 
all spins in the sample: 



5(5 + 1) 



(2) 



The first moment sum rul(f3 provides a direct link be- 
tween raw data and interaction terms in the spin Hamil- 
tonian. 



/C30 
ujS''"'{Q,uj)du 
-OO 



(3) 



1 1 



^Jd(Sr-Sr+d)(l-COsQ.d) (4) 



Here {d} is the set of all bond vectors connecting a spin 
to its neighbors, the index {r} runs over all N spins. The 
Hamiltonian is assumed to take the form 



(5) 



For a powder sample, the magnetic component of 
/(Q, TiLo) is related to the sphericalbi. averaged dynamic 
spin correlation function S{Q,uj) byHd 



(6) 
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FIG. 12. Dispersion of magnetic excitations in CuHpCl de- 
rived from the fits shown in Figs. ^, ^ and 0. Open cir- 
cles represent the lower energy mode, solid triangles repre- 
sent the higher energy mode. Solid lines are phenomenolog- 
ical fits to a dispersion relation satisfying Bloch's theorem: 
huj{Q) = Ao + Ai cos 2Tvh + A2 cos 2nl + A3 cos 27r(/i -I- /)-|- 
J44Cos27r(/l — /). Assuming no mode crossing, the constants 
are 1.00(2) meV, -0.04(2) meV, -0.02(2) meV, -0.07(2) meV, 
and -0.02(2) meV for the lower mode and 1.21(2) meV, 
-0.03(2) meV, 0.00(2) meV, 0.01(2) meV, and -0.03(2) meV 
for the upper mode respectively. 



1.50 




0.5 

C (rlu) 

FIG. 13. Dispersion of magnetic excitations in CuHpCl de- 
rived from the data shown in Figs. ^ and |l^. Open circles 
represent the lower energy mode, solid triangles represent the 
higher energy mode, and the open square point at (101) repre- 
sents the energy of the dominant mode at (100) translated by 
T — (001). Solid lines are phenomenological fits as described 
in the caption to Fig |l^. 



where TZqu> i^Qr^^) is the normahzed instrumentaLres- 
olution functionj£j f{Q) is the magnetic form factorEJ for 
Cu^+j andjg is the average (/-factor, which is g = 2.083 for 
CuHpClEa Carrying out the integration of Eq. || for the 
band of intensity in Fig. ^(a) yields 0.7(1). The proxim- 
ity of the value to S'(S' = 3/4 indicates that this band 
accounts for most of the magnetic scattering from CuH- 
pCl. For a powder sample, the first moment h{uj)Q of the 
measured quantity is related to the spherical 

average of Eq. ^, and is given by 



h{iii)Q = 



ujI„i{Q, 'huj)dhj 



|lf/(Q), ^ 



^ ^ Jd(Sr • Sr+d) ( 1 — ' ^ 



r.d 



Qd 



where d — |d|. The first moment fi{Co)Q computed from 
the data in Fig. ^(a) is shown in Fig. ^ In strongly 
dimerized systems, the term arisiug from the intra-dimer 
bond dominates the first moment In the strongly dimer- 
ized two-leg ladder model of CuHpCl, these are the intra- 
molecular rung bords labeled by 1 in Table 1 and de- 
picted in Fig. |^(a).t3 As shown by the dashed line in 
Fig. 15, the Q-dependence arising from inserting di in 
Eq. 0is manifestly inconsistent with the data. Specifi- 
cally, the di — 3.376 A bond yields a maximum at higher 
Q than is observed in the experiment. Thus, a longer 
bond (or bonds) must give the dominant contribution to 
h{Cj)Q. Fitting the data in Fig. ^ to Eq. ^ with a sin- 
gle, variable bond length da yields da — 7.5(2) A, and 
the dashed-dotted line shown in Fig. |l5[ However, this 
model is still inconsistent with the data. The salient dis- 
crepancy is that the ratio of the low-Q peak intensity to 
that of the high Q plateau is about twice larger in the 
data than in the model. Such a large ratio can only be 
achieved by combining terms of varying sign in Eq. |^. 
With appropriate spin spacings, d, such terms can inter- 
fere destructively in the high Q plateau while the low-Q 
regime is dominated by contributions from long bonds. 
Indeed, the solid line in Fig. |l^, which we shall describe 
in greater detail below, corresponds to a model with both 
positive and negative terms in Eq. |^. According to this 
equation, the magnitude of the high-Q plateau is directly 
proportional to the shift of the ground state energy be- 
low zero while the peak height measures the strength of 
individual spin pair contributions to the ground state en- 
ergy. Hence, there is a direct link between a large peak 
to plateau ratio in first moment data, and frustrated in- 
teractions that raise the ground state energy. 

The single crystal data help to distinguish between the 
eight distinct spin pairs with Cu-Cu spacings in the range 
7 to 7.6 A. The first moment ?i(ci))q of the magnetic scat- 
tering intensity, /m(Q, ^f^), from a single crystalline sam- 
ple is given by 



(7) 




FIG. 14. The combined lattice of spin-spin interactions in 
CuHpCl viewed as a projection on the a — c plane with the 
same color coding as for the individual sub-lattices in Fig. ^. 

^if /(O)P^ E '^d(Sr • Sr+d) (1 - COSQ • d) 



r.d 



Here we have neglected any spin space anisotropy, a rea- 
sonable assumption for a spin- 1/2 quantum spin liquid. 
The first moment may be determined from gaussian fits 
to individual spectra as follows 



;i(c:;)Q = ^/(Q),?iLj,(Q), 



(9) 



h{u!)Q=h / ujI„i{Cl,Tujj)dijj 



where /i(Q) is the integrated intensity for mode i at 
wavevector Q and TiLJi{Q) is the corresponding mode en- 
ergy. The wavevector dependence of ?i(cZ')q is illustrated 
in Fig. ^ where the diameter of the circles is proportional 
to h{uj)Q and as conventional plots along symmetry di- 
rections in Figs. ^ and There is an undetermined 
overall scale factor as the single crystal data was not nor- 
malized in this experiment. The first moment is seen to 
be largest for h « 0.5 along the line {hOl — h). Since the 
magnitude of Q in that part of reciprocal space is close 
to the value 0.55 A~^ where the peak in the first moment 
of the powder data occurs, the powder and single crystal 
data are consistent. The strongest modulation in h{Lo)Q 
occurs along the (hOh) direction, with weaker modula- 
tion also visible in other directions. From the form of 
Eq. ^, this implies that the bond vector that contributes 
most strongly to the first moment is parallel to [101]. 

To determine the relative importance of the magnetic 
bonds, we carried out a simultaneous fit of the powder 
data to Eq. ^ and the single-crystal data to Eq. || with 
a single set of values for the bond energies, Jd(So • Sd). 
Because their contributions to the first moment of the 
scattering data cannot be distinguished, we derive only 
an average correlation term for spin pairs labeled with the 
same numerical index in Table |. Such spin pairs would 
be equivalent had the molecular unit possessed centro- 
symmetry. In addition, there is close similarity between 



> 

CD 

O 

A 

i3 

V 



1.0 


' 4^5 ' 1 ' 1 ' 1 


1 


0.8 




_ 
- 


0.6 


- rXr ''""*» 


- 


0.4 


— f '± / X X X -^[It6f%^i:^r3 
// * J?"Q^ 
#• ' "x" 




0.2 


/' 

" A 







- A / 











0.5 1.0 



2.0 2.5 



1.5 

Q(A-^) 

FIG. 15. First hcj moment of the inelastic powder data in 
Fig. ^ versus wavevector transfer (Eq. The region of in- 
tegration is limited to the band of magnetic excitations from 
0.7 to 1.5 meV. The dashed-dotted line is a fit with a single 
dimer bond length of 7.5(2) A. The dashed line is a fit fix- 
ing the dimer bond length to 3.376 A, corresponding to the 
intra molecular Cu-spacing. The solid line is the best com- 
bined fit to the powder and single crystal data. Fit results 
are described in the text and enumerated in Table | 




FIG. 16. First moment of the magnetic excitations in CuH- 
pCl derived from single-crystal data. Curves are fits as de- 
scribed in text. 



the chemical environments along the exchange pathway 
of a and b labeled bonds. These facts lend some support 
to the assumption that the corresponding bond energies 
are similar. 

The results of the fit are given in the last column of 
Table | and as solid lines in Figs. ^ |l^. It is impor- 
tant to note the direct contribution of a spin pair to the 
ground state energy is small when the exchange constant 
and/or the spin correlation function is small. Further- 
more, according to Eq. ^, negative bond energies lower 
the ground state energy while positive terms indicate a 
frustrated spin pair that raises the ground state energy. 



V. DISCUSSION 

The most remarkable result from Table 1, is that the 
intra-molecular bond 1 is in a frustrated configuration. 
The presence of a frustrated spin pair with a short bond 
vector was anticipated based on inspection of the raw 
powder first moment data. It is easily verified that in 
a closed loop of interacting spins with an odd number 
of antiferromagnetic exchange interactions, only an even 
number of spin pairs can be satisfied. Bond 1 is part of 
no fewer than 5 near neighbor bond triangles that are 
frustrated if all interactions are antiferromagnetic. Fig- 
ure 0(a) shows two of these triangles: (1,3,2) and (1,3,10). 
Fig.l^(b) shows bond triangles (1,4,8) and (1,4,7), while 
Fig. ^(c) shows bond triangle (1,5,6). Of these. Table 1 
clearly indicates that the triang les in Fig. |(b) are frus- 
trated. Our analysis indicates that bonds 4 and 8 are 
satisfied at the apparent expense of bond 1. The energy 
derived for bond 7 while negative, is not statistically sig- 
nificant. This is consistent with the expectation from 



section I A that this is a weak exchange interaction. 

The connectivity of the lattice formed by bonds 1, 4, 
and 8 is that of the gecanetrically frustrated Shastry- 
Sutherland model (SSM)Eil albeit with lower symmetry. 
Compared to the SSM, which can be described as a 
square lattice with alternating diagonal bonds on half 
the squares, the lattice shown by solid lines in Fig. ^(b) 
corresponds to a tetragonal lattice formed by bonds 4 
and 8 with the "diagonal" bonds 1 arranged as in the 
SSM. When the square lattice exchange interactions in 
the SSM are less than 70% of the interactions across the 
diagonal, the groundrState of the SSM consists of singlets 
on every cross bondB For stronger inter-dimer interac- 
tions, there is a first order transition to square lattice 
Neel order. Table 1 indicates that the SSM planes in 
CuHpCl have local spin correlations resembling the Neel 
phase of the SSM. The lattice formed by bonds 1, 5, and 
6 also has the topology of the SSM but correlations in 
this plane do not readily map on a known phase of the 
SSM. 

It does seem surprising that spin pair 1 can provide 
a three times larger positive contribution to the ground 
state energy than the negative contributions from bonds 
4 and 8. To determine whether this is plausible, we ex- 
amined a series of frustrated antiferromagnetic spin-1/2 
clusters. A central spin pair, Sa,Sh with antiferromag- 
netic exchange constant J > is surrounded by 2n spins- 
1/2 Si, S2, ...82,1, which interact with both pair members 
with equal antiferromagnetic exchange, J' > 0. The spin 
Hamiltonian is given by 



Ti — Hq + Ti.2n 

Ho = JSa ■ Sb 



(10) 
(11) 
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FIG. 17. First moment of the magnetic excitations in CuH- 
pCl derived from single-crystal data. Curves are fits as de- 
scribed in text. Data from the (hkO) zone were scaled to data 
from the (hOl) zone using (100) as a point of reference. 



This leads us to the suggestion that bonds 5 and 6 
could play a significant role even if the corresponding 
bond energies appear to be small. As is apparent from 
Fig. 14, bonds 5 and 6 close a molecular triangle that 
projects onto the a — c plane. The description therefore 
emerges of a set of antiferromagnetic SSM layers normal 
to a* coupled in a frustrating triangular lattice geometry 
to a second interpenetrating set of frustrated SSM layers 
normal to b*. A small alternation between the 3a and 
3b bonds that couple layers could be an additional factor 
favoring singlet formation on the stronger of these bonds. 

The magnetic ground state energy may be calculated 
from Eq. ^ and the numbers in Table 1 to be -0.36 meV 
per spin. For comparison, if all bonds were satisfied with 
the same magnitude of spin correlations as observed the 
ground state energy would be -1.17 meV per spin. The 
ratio between the actual ground state energy and the 
latter upper bound on the energy in the absence of frus- 
tration is 0.3. For comparison the ground state energy of 
the SSM at the critical point separating the Neel phase 
and the dimer phase is 0.26 times the energy that the 
spin system would have if all bonds could simultaneously 
be engaged in singlet formation.Ea Hence, CuHpCl is at 
least as frustrated as the SSM at its critical point. 



(12) 



i = l 



We examined clusters with n ranging from 1 to 3. For 
each cluster we diagonalized the spin Hamiltonian and 
determined the ground state frustration index 



/ 



-4n 



< 0|Ho|0 > 

'<0\H2n\0> 



(13) 



as a function of a: = J'/ J in the range where the central 
spin pair is frustrated. We obtained maximum frustra- 
tion indexes fmax = 1,2,3 for n = 1,2,3 respectively. 
The results show that a level of frustration similar to that 
observed for the intra-molecular spin pair in CuHpCl is 
possible even for very small spin clusters. 

While the lattice of Fig. ^(b) provides an explanation 
for the frustration of bond 1, it does not readily account 
for the singlet ground state of CuHpCl. Local correla- 
tions in these planes is Neel-like so the lattice left on its 
own might be expected to have a gapless spectrum and 
long range order at low temperatures. From Table 1 we 
see that bond 3, which generates the lattice in Fig. ||(a), 
is in an un-frustrated configuration. In fact, this bond 
provides the greatest contribution towards lowering the 
ground state energy. Fig. 14 shows how the lattices of 
Figs. 3(a)-(c) intersect to form a three dimensional lat- 
tice. We see that if the lattice in Fig. ^(b) were in a 
Neel phase, then bond 3 only serves to strengthen such 
Neel order and extend it to three dimensions. It is there- 
fore difficult to see how the sub-lattices from Figs, ^(a) 
and ^(b) alone can account for an isolated singlet ground 
state. 



VI. CONCLUSIONS 

In summary, we have presented inelastic neutron scat- 
tering data from deuterated powder and hydrogenous sin- 
gle crystals of the organometallic spin- 1/2 magnet CuH- 
pCl. Consideration of the excitation spectra, the crystal 
structure, and the wave vector dependence of the first 
moment leads to the conclusion that spin-spin interac- 
tions in this system form a complex three-dimensional 
lattice and not a spin ladder as previously thought. While 
structurally one might expect CuHpCl to fall in the spin- 
dimer class of quantum spin liquids, this appears not to 
be the case. The spin gap in spin dimer systems is a con- 
sequence of the singlet ground state of individual dimers. 
However, in CuHpCl the intra molecular spin interac- 
tion is in a frustrated configuration so this cannot be the 
dominant interaction in the problem. Furthermore, if the 
intra-molecular interaction is antiferromagnetic, as ex- 
pected, then the ground state features triplet molecules 
that in and of themselves cannot explain a gap in the 
excitation spectrum. We are therefore led to the conclu- 
sion that the pronounced gap in the excitation spectrum 
of CuHpCl is a consequence of the frustration inherent 
to this particular three-dimensional network of interac- 
tions. This is a surprising discovery as the symmetry 
of the lattice is low. On the other hand, the structure 
clearly is riddled with triangular units, and the connec- 
tivity between them is relatively low. These ingredients 
are known to be important for suppressing Neel order 
and promoting a spin liquid state. E3 Further progress in 
understanding the magnetism of CuHpCl would bene- 



fit from accurate determination of H/D positions using 
neutron scattering, followed by quantum chemical cal- 
culations of exchange constants. More extensive mea- 
surements of the magnetic excitation spectrum are also 
needed, but these must await progress in crystal growth 
or neutron scattering instrumentation. 

There are many organometallic quantum magnets that 
have been labeled as quasi-one-dimensional, based largely 
on the observation of a spin gap. Our experiments on 
CuHpCl have shown how neutron scattering from sin- 
gle crystals as small as 0.1 g can be used to establish 
the dimensionality and the basic nature of interacting 
spin systems. They also show that insulating magnets 
with a gap in their excitation spectrum may constitute a 
considerably more complex and diverse class of interact- 
ing quantum many body systems than previously antici- 
pated. 
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